Mon. Not. R. Astron. Soc. 000, ITHT91 (2011) Printed 26 September 2011 (MN IMfeX style file v2.2) 



O 

(N 



The impact of helium reionization on the structure of the 
intergalactic medium 

Avery Meiksin 1 *, Eric R. Tittley 1 

1 SUPA], Institute for Astronomy, University of Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, UK 
Accepted . Received ; in original form 

ABSTRACT 

We examine the impact of helium reionization on the structure of the intergalac- 
tic medium (IGM). We model the reionization using a radiative transfer (RT) code 
coupled to the combined gravity hydrodynamics code Enzo. Neutral hydrogen and he- 
lium are initially ionized by a starburst spectrum, which is allowed to gradually evolve 
into a power law spectrum over the redshift interval 3.0 < z < 4.0. The temperature- 
density relation of the gas is found to fan out and flatten following Hen reionization, 
with an inversion for high overdensitics of p/(p) > 5. Peculiar velocities of up to 
10 km s _1 are induced by the increased pressure, with the gas density field distorted 
over large coherent regions by 10-20 per cent., and the dark matter by levels of 1 per 
cent. The photoionization-induced flows may thus distort the matter power spectrum 
at comoving wavenumbers k > 0.5 /iMpc -1 by a few per cent, by z — 2. 

Absorption spectra for Hi and Hen are drawn from the simulations, and absorp- 
tion lines are fit to the spectra. The increased temperature following Hen reionization 
reduces the line centre optical depths, resulting in an enhancement in the fraction of 
very low optical depth pixels and an overall distortion in the pixel flux distribution 
compared with Hen reionization in the optically thin limit. A median Doppler param- 
eter of 35km s _1 is obtained for the Hi absorption systems at z = 3. Dividing into 
subsamples optically thick and optically thin at line centre reveals that the optically 
thick systems undergo only mild evolution while the optically thin systems evolve 
rapidly following Hen reionization. A comparison between Hen and Hi absorption 
features shows a broad distribution in the Hen and Hi column density ratio, peaking 
near the measured value and only slightly narrower than measured. The distribution 
of the ratio of Hen and Hi Doppler parameters peaks midway between the thermally 
broadened and velocity broadened limits. A comparison with approximate simulation 
methods based on either a pseudo-hydrodynamical scheme or a non-RT hydrodynami- 
cal simulation with boosted Hen heating rate shows moderately good agreement in the 
absorption line properties, but not to the precision to which they may be measured, 
and not over the full redshift range for which the high redshift Lya forest is observed. 

Key words: radiative transfer - quasars: absorption lines - quasars: general - cos- 
mology: large-scale structure of Universe - methods: N-body simulations 



1 INTRODUCTION 

Cosmological simulations of the intergalactic medium (IGM) 
in the context of cold dark matter (CDM) theories for struc- 
ture formation have proven extremely successful in repro- 
ducing many of the observed properties of the IGM as mea- 
sured through the Lya forest. The mean Hi transmission 
through the IGM is recovered for an Hi ionizing background 
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matching that of the measured Quasi-Stellar Object (QSO) 
luminosity function for z < 3, although an enhanced radi- 
ation field, likely due to galaxies, is required at higher red- 
shifts. Once normalised to the mean transmission, the Hi 
column density distribution is broadly recovered as well for 
a ACDM m odel, as i s the m easured power spectrum of the 
Hi flux. See iMeiksinl [2009) for a review. 

Despite these successes, some of the statistics still show 
lack of agreement with measured values. Prominent among 
these is the line widths of the Hi absorp tion features, which 
tend to be broader than the predictions (Br yan fc Machacekl 
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l200d : iTheuns et ai]|200Ch . A detailed comparison with mea- 
sured spectra reveals that the disagreement is among the 
optically thin Lya systems, for which an additional contribu- 
tion of 10 — 15 km s _1 , added in quadrature to the Doppler 
parameter, is required, while the widths of the optically thick 
syste ms well match the measured values (jMeiksin et al.l 
l200ll ). One possible explanation is recent heat input at z ^ 4 
into the IGM. The additional required broadening corre- 
sponds to a temperature increase of 17 x 10 3 K. The op- 
tically thin systems are associated with moderate density 
or even underdense structures for which the time to reach 
thermal equilibrium exceeds a Hubble time l|Meiksinlll994 
iMiralda-Escude fe Reesl lT994V The optically thin systems 
may thus retain a memory of recent heat input at z 4. 

The most likely source of heating is Hen reionization. 
The abundance of QSOs suggests they would have ion- 
ized Hen at a redshift 3 < ^j£ elI < 4.5, although pos- 
sibly as_ejudy_as zjj en ~ 5 for sufficiently hard spec- 
tra (|Meiksin|[2005l ). Support for recent Hen reionization 
has been based on measurements of the Hen Lya opti- 
cal depth of the IGM. The measurements span the red- 
shift range 2.2 < z < 3.8 based on Hubble Space Tele- 
scope (HST) observations of Q0302 -003 (z cm ~ 3.3) 



jjakobsen et al .11 19941 ; [Heap et alj feOOO). Hopkins Ultraviolet 
Telescope and Far Ultraviolet Spectroscopic Explorer (FUSE) 
observations of HS 1700+ 6416 (z cm ~ 2.7) (|Davidsen et all 
ll99d : lFechner et al J 120061 ). HST and FUSE observations of 
HE 2 347-4342 (z el? ~ 2.9 ) dReimers et al.ll 19971 ; iKriss et all 
120011 ; ISmette et all |2002| ; IShull et all l20ld), and HST ob- 
serva tions of PKS 1935—692 (z em — 3. 2') jAnderson et alj 
1999[) . SDSS J2346-0016 (z em ~ 3.5) dZheng et al.l|2004 
20081 ). Q1157+3143 (z cm ~ 3. 0) (IReimers et al.Bood ) and 
SDSS J1711+6052 (z em ~ 3.8) jZheng et alj200a ). The Hen 
Lya optical depths rise rapidly from thcII — LOO ± 0.07 av- 
eraged ove r the redshift interv al z = 2.2 — 2.6 to THeii — 4.9 
at z ~ 3.3 (Syph ers et al . 2011). as if entering into the epoch 
of Hen reionization. 

The origin of the UV background, however, is still not 
well established. The high Hen Lya optical depths com- 
pared with the corresponding values found for Hi gener- 
ally requires that the UV metagalactic background be soft 
between the Hen and Hi photoelectric edges, with a ratio 
of Hi to Hen photoioni zation rates of \& > 200 at z ~ 3 
(|Madau fc Meiksiniri 994) , corresponding to a source spectral 
index for /„ ~ v~ as of as > 1.8 after all owing for the filter- 
ing of the radiation through the IGM jMadau fc Meiksinl 
1 1994! ; lHaardt fc Madaul 1 19961 ; iBolton et alj l2006h . Thisls 
consistent with the inferred spectral index of bright high 
redshift (0. 3 < z < 2.3) QS Os after correcting for IGM 
absorption (|Telfer et alj l2002h . but it is inconsistent with 
the much harder spectra fou nd for dimmer QSOs nearby 
(z < 0.67) (|Scott et al.ll2004h , for which minimal IGM cor- 
rections are required. Possibly ionizing radiation from galax- 
ies boosts the Hi photoionization rate, however the contribu- 
tion from QSOs alone approximately matches the required 
rate to recover the measu red Hi Lya optical depth (|Meiksinl 
120051 ; IBolton et al.ll2005l ), leaving little room for additional 
sources. 

The optical depth measurements are moreover found 
to fluctuate over a wide range at a given redshift. This is 
most dramatically illustrated by measurements of the Hen 
Lya forest compared with the corresponding measurements 



of Hi. A wide spread is found for the column density ratio 
r\ — Nneii/Nm, var ying at least over the range 4 < r? < 600 
at 2.3 < z < 2.8 (izheng et al] 12004 IReimers etal] 120051 ; 
iFechner et all l200d : lshull et al] 12010 ) . The higher values 
correspond to Hi to Hen photoionization rate ratios of 
* > 1000. Such fluctuations are expected during the epoch 
of Hen reionization, so that the large Hen Lya optical depth 
values may result from a mean over regions in which most 
of the helium is still in the form of Hen. 

The interpretation is not unique, however. Due to the 
discreteness of the sources of ionizing radiation and atten- 
uation by the IGM, local fluc tuations are expected in the 
UV metagalactic backg round (|Zudll992l ; iFardal et al1ll998l ; 
iMeiksin fc Whitel2003l ). The distribution of fluctuations in 77 
measured in HE 2347—4342 are well matched by a UV back- 
ground dominated by QSOs with a soft spectral index con- 
sistent with the data oflTelfer et all d2002h (|Shull et alj2004 
IFardal et al"1ll998l ; IBolton et al.ll200d), as are th e r\ fluctua- 
tions measured in HS 1700+6416 l|Meiksinll2009l ). Possibly a 
harder spectrum could be accommodated in the presence of 
additional Hi ionizing sources other than QSOs, but, again, 
there is not much room to add more sources without over- 
ionizing the hydrogen. The full range of fluctuations in the 
measured values of THeii in several QSO lines of sight sug- 
gests the epoch of H en reionization may have completed as 
recen tly as z ~ 2.9 l|Furlanetto fc Dixonlfeoid ; IShull et al] 
2010). By contrast, recent attempts to quantify the thermal 
evolution of the IGM from the Hi Lya forest measured in 
QSO spectra indicate a ris e in the temperatu re of gas at 
the mean density for z ~ 4 (|Becker et al][201ll ). suggesting 
Hen reionization was becoming widespread at these times. 
If both interpretations are correct, then Hen reionization 
was a drawn out process extending over the redshift range 
2.9 < z < 4.5. 

The purpose of this paper is to investigate the dynami- 
cal impact of helium reionization on the structure of the in- 
tergalactic medium and its observational signatures. While 
early simulations of the IGM included helium reionization, 
they were done in the optically thin limit. Cosmological sim- 
ulations solving the radiative transfer (RT) equation during 
helium reionization were restricted to the gravitational com- 
ponent, wh ich was used to scale the baryonic componen t 
properties (|Tittlev fc Meiksinl l2007i ; iMcQuinn et~all 12000 ) . 
In this paper, w e implement radiatio n hydrodynamics using 
the algorithm of lBolton et alj (|2004l ) coupled to the gravity- 
hydrodynamics code EnzcQ (specifically v. 1.0.1). The hy- 
drodynamical response of the gas to the boost in heating 
when radiative transfer is included could result in more rapid 
photo-evaporation of the gas in small haloes, allowing reion- 
ization to occur more rapidly, as well as in observational 
signatures on the Lya forest such as increased line broad- 
ening due to outflows from the haloes and a decrease in the 
ratio of gas density to dark matter density. An increase in 
the gas temperature will also alter the ionization fractions of 
hydrogen and helium at a given gas density, producing mod- 
ifications to the column densities of the absorption features, 
the distribution function of pixel fluxes, and the flux power 
spectrum, all basic statistics used to quantify the Lya forest 
and the predictions of cosmological models for its structure. 
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All results are for a flat ACDM universe with the 
cosmological parameters Q, m = 0.24, Qth 2 = 0.022 and 
h = Ho/WO km s _1 = 0.73, representing the total mass 
density, baryon density and Hubble constant, respectively. 
The power spectrum has spectral index n = 0.95, and is 
normalized to <J 8h -i = 0.74. 

This paper is organized as follows. Estimates for the ex- 
pected IGM temperature boost following Hen reionization 
are discussed in the next section. The simulations are de- 
scribed in Sec. 3 and the results presented in Sec. 4. The 
Hi and Hen spectral signatures of Hen reionization are pre- 
sented in Sec. 5. A comparison with approximate simulation 
methods is provided in Sec. 6. The principal conclusions are 
summarised in the final section. 



2 PHOTOIONIZATION HEATING 

A radiation field with its specific energy density locally 
approximated near the Hen photoelectric edge by u v — 
Uhiy '/vL,)~ a will photoionize the Hen at the rate per Hen 
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where vl is the frequency of the Hen Lyman edge and 
the photoionization cross-section is approximated as a ~ 
croiv /vl)~ 3 '. The corresponding heating rate per Hen ion is 

G Hen = / Av^ihu ~hv L )a v 
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The heating rate per ionization, normalized by the ionization 
potential, is then e Hen = G HelI /(/i^zTji elI ) = 1/(2 + a). 
The energy injected into the gas is thus sensitive to the 
shape of the local ionizing spectrum. A hardened radiation 
field with a < can result in a large amount of energy 
deposited per ionization. 

The spectrum is expected to harden within an ioniza- 
tion front because hard photons are less likely to be ab- 
sorbed than soft photons above the photoeletric threshold 
for a given optical depth at the threshold energy. In gen- 
eral, for a column density N the optical depth above the 
photoelectric threshold of hydrogen or singly ionized helium 
is t u — tl(v/vl)~ 3 where tl = aoN. The photoionization 
rate per atom/ion may be expressed in terms of the incident 
radiation density from the source uf, = uf(v/vL)~ as as 
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where y(a,x) = f£ dte t t a 1 is the incomplete gamma 
function. Similarly, the heating rate per atom/ion is 
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Figure 1. Temperature increment for ionizing Hen to He III as 

a function of the optical depth T£ at the Hen Lyman edge. 

Shown for a source spectrum f v oc v~ a s with spectral index 

ag = 1.8 (solid line), 1.5 (long-dashed line), 0.5 (short-dashed 

line), (dot-dashed line) and -0.5 (dotted line). The upper set 

of curves corresponds to reionization going to completion at a 

fixed optical depth, as would occur for radiation filtered through 

a foreground absorption system in ionization equilibrium. The 

temperature boost increases for harder source spectra. For large 

1 /3 

tl, the temperature boost increases like t l . The lower set of 
curves corresponds to the temperature boost for a time varying 
optical depth as an ionization front sweeps through a slab of ma- 
terial with an initial optical depth r^(0) given by the value 
indicated. Very high temperatures are reached for hard spectra. 
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The resulting heating rate per ionization of Hen will 
produce a temperature increment 
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where Ti is the initial gas temperature and £j£ elI = 
G/(hviX) for Hen. The increment corresponds to the com- 
plete ionization of Hen to He in at fixed tl, as would occur in 
the presence of an intervening Lyman limit system, already 
in ionization equilibrium, lying between the source and the 
region being photoionized if the Lyman limit system domi- 
nated the optical depth at the Lyman edge. The radiation 
field reaching the region being photoionized would be hard- 
ened by the absorption within the intervening system. The 
temperature increment is shown as a function of the opti- 
cal depth tl at the Hen Lyman edge for a range of source 
spectral indices in Fig. [T] (upper set of c urves). A helium 
mass fraction of Y = 0.248 was adopted ( Stei gmanl I2007T ) . 
A system with overdensity p/{p) and comoving thickness L 
would have an optical depth at the Hen Lyman edge of 
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The radiation hardening through overdense systems with 
P/{p) ~ 200 and comoving sizes of 0.05 — 0.25 kpc would 
produce temperature boosts in the "shadows" they cast in 
less dense, more quickly photoionized structures, of up to 
AT ~ 40000 K, and even reaching 10 5 K in the shad- 
ows of intervening systems with TL,Heii > 100. Radiative 
recombination and excitation atomic line cooling will sub- 
sequently lower the temperature of the gas once the ion- 
ization front has passed, although the equilibrium temper- 
ature may never be achieved in rarefied und erdense regions 
because of the long require d time scales l|Meiksinl 1 19941 : 
iMiralda-Escude fc Reesll 19941 ). 

Alternatively, the heating produced by the passage of 
an ionization front through a slab of gas in which the optical 
depth decreases as the ionization front sweeps through may 
be approximated by equating the incident rate of ionizing 
photons with the rate at which the gas in the slab is pho- 
toionized. For an initial optical depth tl(0) at the Lyman 
edge, the optical depth at the edge will evolve according to 
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where k = 1/as for as > 0, k = log(^ max /^L) for as = 
and k = [(^ m a X /^L)~ QS - l]/(— as) for a s < 0, where ^ max 
is the maximum frequency for which photoionizing photons 
are produced. Using Eqs. Q and Q for the photoionization 
and heating rates, the number density n of the species un- 
dergoing ionization and the thermal energy input per unit 
volume u, as a function of tl, are governed by 
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The resulting post-ionization temperature is shown in Fig. [T] 
(lower set of curves). Temperature boosts of ~ 1000 — 
10000 K are typical, although hard spectra sources (as < 0) 
can produce substantially higher temperatures. (Somewhat 
arbitrarily, v ma x/vL = 4 was assumed for as < 0.) Again, ra- 
diative recombination and excitation line cooling will even- 
tually decrease the gas temperature. 

Temperature boosts following Hen reionization typi- 
cally ranging between AT ~ 5000 - 20000 K at z = 3, 
but reaching as high as AT ~ 30000 - 40000 K directly af- 
te r reionization, were found in the reionization simulations 
of iTittlev fe Meiksin! (|2007T ) for a source with spectral index 
as = 0.5, turning on after a starburst-like spectrum pho- 
toionizes the hydrogen and neutral helium to singly-ionized 
helium. For a pure power-law source spectrum (as = 0.5) 
and for a mini-quasar spectrum, lower overall temperatures 
of 5000 - 20000 K were obtained at 3 < z < 4. Higher 
temperatures were generally found in the denser structures. 
Te mperature boost s of A T ~ 10000 - 20000 K were found 
by iMcQuinn et al l (120091 ) for a QSO spectral distribution 
with a mean spectral index of as = 1.2, and boosts of up to 
AT ~ 30000 K were found for harder source spectra with a 
mean spectral index as = 0.6. 



In addition to cooling by atomic processes, the gas will 
undergo adiabatic cooling losses as it escapes from gravita- 
tional potential wells that were able to bind the gas when 
cooler. The expan sion velocities will also broaden the ab- 
sorption features |Meiksinl 1 19941 ). Previous Hen reioniza- 
tion simulations, which used the dark matter distribution 
to model the baryon density, were not able to take these 
effects into account. In this paper, we take the hydrody- 
namical response of the gas into account by coupling our 
radiative transfer scheme to a gravity-hydrodynamics code. 



3 SIMULATIONS 

We have modified the combined hydrodynamics iV-body 
gravity code Enzo by replacing its ionization and atomic 
cooling modules by ones employing the methods used by 
ITittlev fc Meiksin! (|2007l ). In summary, a probabilistic ra- 
diative transfer method is used to compute the reionization 
of both hydrogen and helium. Atomic cooling from the ra- 
diative recombination of the hydrogen and the two helium 
states is included, along with collisional excitation of neu- 
tral hydrogen and Compton cooling. Cooling due to the col- 
lisional ionization of hydrogen and the collisional excitation 
and ionization of helium are negligible for the temperatures 
encountered in the simulations. 

For the hydrodynamical processes at the resolution of 
cosmological IGM simulations, the dynamical time scales are 
much longer than the radiative time scales in the vicinity 
of the ionization front. It would be prohibitively expensive 
computationally and unnecessary to force the hydrodynam- 
ics to evolve at the same time scale as the radiative trans- 
fer. Hence, in our implementation the radiative transfer time 
scale is allowed to drop below that of the hydrodynamics. 
The energy equation is solved through operator splitting by 
computing the ionization state of the gas and the amount of 
•(photoionization heating and atomic cooling between the hy- 
drodynamical time steps. During the ionization timesteps, 
all parameters of the gas are taken to be those at the end of 
the immediately preceding hydrodynamical timestep with 
the exception of the fluid density and temperature. The 
density is varied linearly from the previous hydrodynamical 
timestep to that of the current ionization timestep while the 
radiative transfer processes along with the gas temperature 
are evolved on their own, shorter, adaptive time scale. 

The simulations were carried out in a box 25fe _1 Mpc 
(comoving) on a side. This corresponds to the order of the 
expected size of a Hen ionization front at the helium reion- 
ization epoch, when the characteristic diameter of the He in 
zone produced by an individual QSO is comparable to the 
mean distance betwe en QSOs, which i s typically ~ 5 Mpc 
(comoving) at z ~ 3 | Wolf et al.ll2003l ; lMeiksTnll2005l ) . Then 
Enzo simulations were performed with 512 3 gravitating par- 
ticles and 256 3 grid zones for the fluid component, corre- 
sponding to a resolution of 98 kpc (comoving), sufficient to 
resolve the Jeans length. 

The incident spectru m corresponds to a h ybrid 
starburst-QSO model as in ITittlev fc Meiksin! (|2007h : an 
initial starburst spectrum is turned on at z = 8 and trans- 



formed into a power-law spectrum /„ oc v 



turning on 



the power-law spectrum initially at z — 4 and ramping it 
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Figure 2. The heating rate per Hell ionization at z = 3.3 for a 
Hen ionization front sweeping across the box from top to bottom. 

up to its full value at z — 3, while the starburst spectrum is 
ramped down during the same interval. 



4 RESULTS 

4.1 Physical impact 

The intensity of the starburst spectrum was adjusted to en- 
sure the hydrogen ionization front swept across the simu- 
lation volume by z ~ 6, although some shadowed regions 
persist along the edge of the box until z ~ 5.5. The in- 
tensity of the power-law spectrum was adjusted so that the 
Hen-ionization front swept across the box by z ~ 3.2, again 
with some shadowing persisting to z 3. This corresponds 
to a hydrogen ionization rate from the QSO at full intensity 
half that of the starburst at its full intensity. 

The heating rate per ionization of Hen as the Hen ion- 
ization front passes through the box at z = 3.3 is shown in 
Fig. [2] The filtering effect of the IGM in hardening the ion- 
izing radiation field is clearly discernible in the pre-ionized 
gas ahead of the ionization front. The heating rate per ion- 
ization is substantially enhanced, reaching values of up to 
100 — 400 eV. By comparison, once the front has passed, the 
energy per Hen ionization is ~ 22 eV. 

The resulting temperature structure is shown in Fig. [3] 
as the Hen I-front passes through the box at z = 3.3 and 
at z — 2.0, after Hen ionization has completed. The effect 
of shadowing produced by density inhomogeneities in the 
IGM are apparent as the advance of the ionization front is 
delayed beyond the densest regions. This may best be seen 
by comparing the temperature map at z = 3.3 in Fig.[3]with 
the density map in Fig. [4] (The principal density structures 
change little between z = 3.3 and z = 2.0 in comoving 
coordinates.) Large temperature boosts of 3 — 4 x 10 4 K 
occur downstream from dense clumps, consistent with the 
striations of high heating rates per ionization in Fig. [5] as 
expected for photoionization by a radiation field hardened 
by intervening optically thick clumps with optical depths of 
several tens, as discussed in Sec. [5] 

By z = 2.0, most of the high temperature striations 
have been replaced by gas temperatures that trace the mat- 
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Figure 5. Distribution of temperatures at z = 3 following reion- 
ization including radiative transfer (solid line) vs reionization in 
the optically thin limit (dashed line). 

ter distribution now that Hen reionization has completed. 
Pronounced hot regions, however, persist, particularly in 
underdense regions unable to achieve thermal equilibrium. 
These lay behind dense structures while the gas was reion- 
ized. The filtering of the radiation field by the dense struc- 
tures hardened the spectrum of the incident radiation pass- 
ing downstream to the underdense gas, as shown in Fig. [2] 
producing the high temperatures. This effect is especially 
prevalent in the centre-right portion of the figure. 



4.2 Differential impact 

Including radiative transfer in the ionization of Hen sub- 
stantially boosts the temperature of the gas, increasing it 
typically by 15 - 20 x 10 3 K over the optically thin reion- 
ization limit, as shown in Fig. [5] The temperature-density 
relation is shown in Fig. [6] The plateau with a gas tem- 
perature of ~ 20 — 30 x 10 3 K at th e mean baryon den - 
sity agrees well with the estimate of ISchave et alj (|2000h 
of T ~ 20 - 30 x 10 3 K at z = 3, with only a weakly 
increasing temperature with overdensity, as inferred from 
the Doppler widths measured in high resolution spectra of 
the Lya forest. The temperature found in the simulation, 
however, is much larger than the temperature at the mean 
density of To ~ 10 — 15 x 10 3 determined from the flux cur- 
vatur e method applied to high resolution data dBecker et alj 
l201ll ). The discrepancy may be due to an inadequate sim- 
ulation model for calibrating the statistic used to estimate 
the gas temperature. Alternatively, it may indicate that ei- 
ther Hen was reionized too late in the simulation here, or 
that the Hen-ionizing source was too hard. A temperature 
inversion occurs for pj (p) > 5, as is expected for thermal bal- 
ance between photoionization heating and atomic line and 
radiative recombination cooling in dense structures. Only 
such high density gas is able to achieve thermal balance; 
the time scal e to achieve equilibrium is too long in un - 
derdense gas (lMeiksinlll994l ; iMiralda-Escude fc Reed 1 19941 ) ■ 
A similar trend between temperature and density is found 
for a pure dark matter simulation with radiative transfer 
l|Tittlev fc Meiksinll2007h . 

Compared with optically thin reionization, the excess 
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Figure 3. The post-ionization temperature for a Hell ionization front sweeping across the box from top to bottom, at z = 3.3 (left 
panel) and z = 2.0 (right panel). The box side is 25/i~ 1 Mpc (comoving). 



23.5 




23.5 




Figure 4. The gas distribution at z = 2.0 without Hell rcionization (left panel) and with (right panel). Gas is driven out of small haloes 
by the Hell rcionization. The box side is 2hh~ 1 Mpc (comoving). 



temperature boost allowing for radiative transfer occurs pri- 
marily for underdense g shown in Fig. [S] (right panel). 
Most of the gas with p/(p) > 3 experiences very little boost 
at all, as it is able to quickly recover thermal balance. 

A map of the difference in the temperatures is shown in 
Fig. Boosts of AT > 15 x 10 3 K are visible downstream 
of very dense structures which act to filter and harden the 
radiation field passing through them. The densest struc- 
tures themselves show little change in temperature, as the 
temperature rapidly reaches thermal balance. Small over- 
dense patches, however, show marked temperature deficits, 
as much as AT < — 10 x 10 3 K. These regions account for the 



temperature reductions in Fig. [6] The temperature decreases 
arise from complete shadowing of the Hen- ionizing photons 
by an intervening optically thick system at the Hen photo- 
electric edge. While for optically thin reionization, all the 
gas everywhere is exposed to the photoionizing background, 
allowing for radiative transfer will produce shadowed regions 
in the presence of sufficiently dense clumps of gas. Only af- 
ter Hen reionization completes will the regions be reheated 
by overlapping ionization fronts. 

The increase in the peculiar velocity compared with the 
computation in the optically thin approximation is shown in 
Fig. [8] (right panel). The differences are small, typically less 
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Figure 6. (Left panel) The temperature as a function of gas density at z = 3.0 following Hen reionization including radiative transfer. 
(Right panel) The boost in temperature as a function of gas density at z = 3.0 including radiative transfer vs reionization in the optically 
thin limit. The contour levels are probability density levels (per dT-dlog 10 p/(p)) stepped by 0.5 dex starting at 10 — 7 . 



25 km/s EnzoRT - Enzo z=3.0 25 km/s EnzoRT - EnzoRT-OT z=3.0 
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Figure 8. Differential peculiar velocity map in a kinematically complex region at z = 3.0. (The indicated box length scale is comoving.) 
Left panel: Map showing the change in the projected peculiar velocity of the gas with Hen reionization including radiative transfer 
from the case without Hen reionization. Right panel: Map showing the change in the projected peculiar velocity with Hell reionization 
including radiative transfer from the case with Hell reionization in the optically thin approximation. 



than 1km s" 1 as shown in Fig. [9] (dashed curve), except 
near the most overdense structures, or behind overdense re- 
gions which shadow the incident radiation field and delay 
full helium reionization. The small magnitudes are in con- 
trast to the peculiar velocity boosts induced by reionization 
compared with the non-reionization case, shown in the left 
panel. Reionization itself induces sizeable peculiar veloci- 
ties, typically 0.1 < Av < 10 km s _1 , throughout most of 
the overdense g shown in Figs. [S] (left panel) and 

(solid curve). 

A comparison of the two panels in Fig. [4] shows that 
some of the dense gaseous clumps present in the run with- 



out Hen reionization disperse when Hen reionization is in- 
cluded, particularly in the complex filamentary regions in 
the lower half of the figure. The corresponding boosts in 
the outflows following Hen reionization are shown in the 
left panel of Fig. [5] The differential impact of Hen reion- 
ization on the gas density of the IGM is shown in Fig. QIJ] 
(left panel). Substantial changes are found, with a ~ 30 per 
cent, decrease in underdense regions as the increased pres- 
sure gradients resulting from the additional heat input drives 
gas out of the shallower potential wells. Enhancements are 
produced as well, with the density increasing by as much 
as 50 per cent, in dense regions when helium reionization is 
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Figure 10. Maps showing the differential effect of Hell reionization on the gas density at z = 3. Left panel: Map showing the percentage 
change in the gas density with Hell reionization including radiative transfer from the case without Hell reionization. Right panel: Map 
showing the percentage change in the gas density with Hell reionization including radiative transfer from the case with Hen reionization 
in the optically thin approximation. (Note change in scale.) The box side is 25h~ Mpc (comoving). 




Figure 11. Maps showing the differential effect of Hen reionization on the dark matter density at z = 3 (left panel) and 2 = 2 (right 
panel). The box side is 25h _1 Mpc (comoving). 



included. These regions tend to trace the filamentary struc- 
ture of the gas distribution, suggesting they are too hot to 
accrete as readily on to nearby haloes. 

Allowing for Hen reionization, but in the optically thin 
limit, produces a similar gas density field to that with ra- 
diative transfer, as shown in Fig. [10] (right panel). Large 
differences, however, occur in the more complex density re- 
gions, with the radiative transfer case generally resulting in 
more completely evacuated underdense regions. 

The gas motion following Hen reionization also pro- 
duces a change in the dark matter distribution, although 
at a much smaller level. As gas is driven out of filamentary 



structures, the reduced gravitational potential will result in 
the dark matter readjusting to a somewhat lower density, 
as shown in Fig. 1111 The differences are small at z — 3, on 
the order of 1 percent or less. By z = 2, however, the dif- 
ferences have grown to as high as 3 per cent., with coherent 
structures several comoving megaparsecs in scale showing 
underdensities. As shown in Fig. 1121 most of the differences 
are accounted for in the optically thin limit, but not en- 
tirely, especially in complex density regions such as in the 
lower middle region of the figure. Large-scale simulations 
with multiple QSO sources are necessary to estimate the 
impact on the matter power spectrum P(k), however the 
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Figure 12. Maps showing the differential effect of Hell reionization on the dark matter density allowing for full radiative transfer 
compared with Hen reionization in the optically thin approximation. Shown at z = 3 (left panel) and z = 2 (right panel). The box side 
is 25h~ 1 Mpc (comoving). 
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Figure 7. Map showing the change in the temperature including 
radiative transfer and without, at z = 3.0, following Hell reion- 
ization. The box side is 25 h~ 1 Mpc (comoving). 



results here suggest reionization may perturb P(k) by a few 
per cent, on comoving scales k > 0.5/iMpc -1 by z = 2. 



5 SPECTRAL SIGNATURES 

To assess the impact of radiative transfer on the spectral sig- 
natures of the Lya forest, spectra against a ficticious back- 
ground source are drawn from the Enzo simulations with and 
without radiative transfer. The s pectra are generated follow- 
ing the procedure described in iTittlev fc Meiksinl (|2007h . 
They are constructed at an angle relative to the simula- 
tion axes in order to create a long non-repeating spectrum 
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Figure 9. Distribution of peculiar velocity difference at z = 3 
for reionization including radiative transfer vs the case with no 
reionization (solid line) and reionization in the optically thin limit 
(dashed line). Also shown is the total peculiar velocity distribu- 
tion for the reionization simulation including radiative transfer 
(dot-dashed line). 



Table 1. Effective IGM optical depths for Hi and Hen Lya. 
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that wraps around a slice through the simulation volume. 
The velocity pixel widths are 5km s _1 . Both Hi Lya and 
Hen Lya spectra are made. The spectra are normalised to 
a mean transmission of (exp(— t)) = cxp(— T e s) with the 
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Figure 13. (Top left panel) The pixel flux distribution function 
at z = 3 following Hen reionization for Enzo with radiative trans- 
fer (black solid line) and in the optically thin limit (red solid line). 
Also shown is the recovery of the flux distribution using only the 
fit absorption features for the run with radiative transfer (blue 
short-dashed line) and without (cyan short-dashed line). (Bottom 
left panel) The corresponding cumulative flux distribution. (Top 
right panel) The distribution of line centre optical depths for the 
fit lines for reionization including radiative transfer (black solid 
line) and without (red long-dashed line). (Bottom right panel) 
The corresponding cumulative distribution. 



values giv en in Table [J for Hi ([Faucher-Giguere et al.|[2008l ) 
and Hen (|Zheng et alj|2004l ; IShull et alj|2010h. Absorptio n 
lines are fit to the spectra using AUTO VP ilbave et al1ll997l) , 
modified to recover from fatal errors (jMeiksin et al.l 12001). 
At z = 3.0, the Hen Lya transmission is very patchy and 
found too small to readily fit absorption line features. 



5.1 Hi absorption signature 

The pixel flux distribution at z = 3.0 is shown following 
Hen reionization both with and without radiative transfer in 
Fig. [TH] (left panels). Allowing for radiative transfer produces 
a smaller fraction of low optical depth pixels (0.1 < r < 0.6), 
and a larger fraction of very low values (r < 0.1), compared 
with the optically thin approximation (right panels). 

Also shown are the flux distributions recovered by fit- 
ting absorption lines to the spectra, and regenerating the 
flux distribution from the absorption lines alone. These are 
found to accurately recover the original flux distribution, 
demonstrating the absence of residual absorption that may 
not be accounted for by absorption lines. 

The distribution of line centre optical depths is provided 
in the right panels of Fig. 1131 showing a higher fraction of 
lines with line centre optical depth to < 0.1 for the case with 
radiative transfer compared with the optically thin limit. 
The higher proportion of very low fluctuations when radia- 
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Figure 14. (Top left panel) The Hi column density distribution 
at z = 3 following Hell reionization for Enzo with radiative trans- 
fer (solid line; black) and in the optically thin limit (dashed line; 
blue). (Top right panel) The corresponding Doppler parameter 
distribution. (Lower left panel) The Doppler parameter distribu- 
tion for absorption lines with line centre optical depth tq < 1. 
(Lower right panel) The Doppler parameter distribution for ab- 
sorption lines with tq > 1. 



tive transfer is included may appear counter-intuitive since 
the temperatures, and so thermal pressures, are greater. The 
increased Jeans length may have been expected to smooth 
away small fluctuations in the spectra. Consideration of the 
dependence of the optical depth on the thermal width pro- 
vides an explanation. Since the line centre optical depth 
varies as to cx Nm/b, for a given column density, an in- 
creased Doppler parameter due to heating the gas to a higher 
temperature reduces the line centre optical depth. The in- 
crease in temperature will also decrease the recombination 
rate at a given gas density, and so reduce the Hi column 
density. 

The column density and Doppler parameter distribu- 
tions for the absorption line fits to the spectra are shown 
in Fig. 1141 The column density distributions for the simula- 
tions with and without radiative transfer are very similar, 
with the distribution for the RT case somewhat broader. By 
contrast, the Doppler parameter distribution for the RT case 
shows a higher cutoff at the low end compared with the op- 
tically thin case. Dividing the lines into those optically thin 
and optically thick at line centre reveals that a large contri- 
bution to the offset originates from the optically thin lines. 
The median Doppler parameter of the optically thin absorp- 
tion lines for the RT case is 6 me d — 45 km s~ , while for 
the optically thin case, 6 me d — 41km s _1 . For the optically 
thick absorbers, 6 mc( j — 26 km s _1 for the RT case, while for 
the optically thin reionization Ca,Se &mcd ~ 24 km s 1 . Since 
the optically thin absorbers arise in moderate to low over- 
density structures, the gas in the optically thin systems is 
too rarefied to establish thermal balance between radiative 
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Figure 15. (Top panel) The Hi Doppler parameter cumulative 
distribution for Enzo with radiative transfer, shown at z = 2 (solid 
line; black), 2.5 (long-dashed line; blue), 3.0 (short-dashed line; 
green) and 3.5 (dotted line; red). (Bottom panel) As in the top 
panel, but with the lines divided between those with line centre 
optical depth to < 1 (rightmost curves) and those with to > 1 
(leftmost curves). 



Figure 16. (Top panel) The ratios of the effective optical depths 
n.Ly^/n.Lyo (solid points) and T; jLy7 /r iiLya (solid squares) for 
Hi obtained from the Enzo simulation with radiative transfer, 
as a function of redshift z. Also shown is the predicted ratio 
T l,Lyy/ T l,hya (open squares) based on a line-blanketing model 
(see text). (Bottom panel) The characteristic Hi T^ ya optical 
depth inferred from the values of r; Lys/^Lya from the Enzo 
simulation with radiative transfer. 



heating and at omic cooling. As a consequence, the systems 
are overhe ated (Miralda -Escude fc Reeslll994l ; lMeiksinlfl99l 
iMeiksin et al.ll200ll ). 

A similar tre nd was found for the optically thin reion- 
ization models of IMeiksin et al.l ([2001) compared with the 
data: the median Doppler parameter of the measured opti- 
cally thin lines were broader than the lines drawn from the 
best-fitting ACDM model simulation. The median Doppler 
parameter fo r lines with tq < 1 mea sured in the spectrum of 
Q 1937-1009 jBurles fc Tvtledl 19971 ). over the redshift range 
3.1 < z < 3.7 was found to be 6 me d — 29 km s _1 , while the 
best-fitting model gave & me d — 20 km s _1 . The measured 
median value for the systems with ro < 1 is substantially 
smaller than found here for the RT simulation, suggesting 
the Hen reionizat i on wa s too recent in the si mulation (cf. 
iTittlev fc Meiksinl (|2007l ); IMeiksin et all l|2010l) ). or that the 
input spectrum too hard on average. By contrast, the data 
for the ro > 1 systems gave 6 mc( j — 26 km s _1 , in agreement 
with the RT simulation result found here. A fully fair com- 
parison with the observations, however, would require repli- 
cating the redshift range, resolution and noise properties of 
the measured spectra, as well as varying the luminosities, 
spectral shapes and timings of the sources. Such a compar- 
ison would require a large suite of simulations, which is not 
the purpose of this paper. 

The evolution of the Doppler parameter cumulative dis- 
tribution is shown in Fig. 1151 The distributions show little 
evolution over 2.5 < z < 3.5, but a marked decrease in 
the Doppler parameters by z — 2.0 as the IGM adiabati- 
cally cools following Hen reionization. Dividing the absorp- 



tion lines into those optically thick and thin at line centre 
(bottom panel), reveals that the optically thick systems un- 
dergo continual evolution towards smaller values with cos- 
mic time, with the median Doppler parameter evolving from 
fe mc d ~ 33 km s _1 at z — 3.5 to 19 km s _1 by z = 2.0 as 
the gas establishes thermal balance following Hen reioniza- 
tion. By contrast, the systems optically thin at line centre 
show almost no evolution between z — 3.5 and 3.0, with 
heavily broadened lines following Hen reionization having a 
median Doppler parameter fo mo d — 44 km s _1 , diminishing 
slightly to 39 km s^ 1 at z = 2.5, and then substantially to 
27km s" 1 by z = 2.0. 

In addition to the scattering of Lya photons, the higher 
order Lyman resonance line photons will be scattered as 
well. Identifying individual absorption features in spectra 
is difficult except for the strongest lines, especially at high 
redshifts. The effective optical depths resulting from line 
blanketing, however, may be inferred. The observed effective 
optical depth at redshift z due to systems with restframe 
equivalent widths w (in wavelength units) is given by 

1 + 2 r d 2 N 

n,-Lyn = -: / aw W, (10) 

Ain J OZOW 

where Ai n is the wavelength of the Lyman resonance line 
transition to the n th principal quantum level, and where 
the equivalent width is related to the line centre optical 
depth ro, in by w = 2(b/c)\i n F(TQ,i n ). . where F(t) is a func- 
tion dependent on the Voigt profile ([Mciksin 200^). If the 
line distribution is separable into the form d 2 N/dzdw = 
(dN/dz)f(w/wl n ), where the equivalent width distribution 
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Figure 17. (Top panel) The distributions of the ratio of the pixel 
optical depths ^Hen/Tii ^ n the spectra following Hell reioniza- 
tion for Enzo with radiative transfer at z = 2.0 (solid line; black), 
2.5 (long-dashed line; blue) and 3.0 (short-dashed line; red). (Bot- 
tom panel) The corresponding cumulative distributions. 



f(w/wl„) may be characterised by a single parameter w* n , 
such as the scale height for an exponential distribution 
f(w/wl n ) = exp(—w/w* n ), then the optical depths for tran- 
sitions to levels n and n' are related by 



Tl.n' 
Tl,n 



F 



A, 



- r 0,lr 



(11) 



where r * ln is given by w* n = 2(b/c)X ln F (r * ln ). 

The ratios r ;iLyi g/r( iLyci and r ;iLy7 /r iiL ya from the Enzo 
simulation with RT are shown in the upper panel of Fig. 1161 
A rise in both ratios occurs at redshifts z < 4 once Hen 
reionization begins. The inferred values of the characteristic 
Lya optical depth r^ ya are shown in the bottom panel. The 
corresponding predicted values of r(.L y7 /T; i Lyci (upper panel) 
lie close to, but somewhat below, the values measured from 
the simulation, suggesting the single parameter separable 
model for the line distribution provides a good description 
of the line blanketing, but becomes less accurate at large 
values of r£ ya . 

5.2 Hen absorption signature 

For a homogeneous radiation field, the ratio of the Hen to 
Hi optical depths per pixel is 



1 AhcII 

4 



&HI 



(12) 



THcII 

THI 4 jVm &HcII ' 

where Nm and JVHeii are the Hi and Hen column densities, 
respectively, and 6hi and &Heii the corresponding Doppler 
parameters. For pure thermal broadening, &h c ii = bm/2, 
while for velocity- broadened lines, i>Hcii = bm- For meta- 
galactic Hi and Hen photoionization rates Thi and Timi, 



the column density ratio, for a hydrogen to helium number 
density ratio of 12.9, is 



Tni "-He QHoIII 
ThcII "h Qhi 

r H i 7.107 - (1/2) logT + 0.00547T 1/3 



0.31 



0.42 



FhcII 6.414 

r H i 



(1/2) \ og T + 0.00868T 1 / 3 



(13) 



where the last form is accurate to 5 per cent, for 1 4 < T < 
5 x 10 4 K |Madau fc Meiksin|[l99i lMeiksinll2009l '). The op- 
tical depth ratio is thus a measure of the fluctuations in the 
ratio of the radiation field as well as the velocity-broadening 
of the absorption lines. 

The distribution of the optical depth ratio per pixel is 
shown in Fig. [IT] for 2 < z < 3. The distributions show 
a wide range of values, with the peak declining from 105 
at z = 3 to 75 at z = 2. The FWHM of the distribu- 
tions is about 30, changing little with redshift. The varia- 
tions arise entirely from the radiative transfer of the incident 
radiation through an inhomogeneous medium following he- 
lium reionization. Allowing for the Poisson fluctuations from 
the finite number of sources within an attenuation volume 
for H e ii- ionizing photon s will further enhance the fluctua- 
tions dFardal et al.lll99Sl; iBolton et alj|2006l ; lMeiksirJ[2009l ; 
iFurlanetto fc Dixonll2oTfj| ). 

Absorption lines were fit in the Hen spectrum at z = 
2.5, and matched to the Hi features, requiring a line-centre 
velocity difference smaller than 10 km s _1 for a successful 
match. The relation between the matching line centre optical 
depths is shown in Fig. 1181 While the deblending of absorp- 
tion features will introduce scatter into the relation, a clear 
correlation is apparent in the upper left panel corresponding 
to To,Heii/To,Hi ~ 100. The ratio is found to decrease towards 
increasing to, hi (lower left panel). This arises as the Hi fea- 
tures become sat urated. A similar tren d was found for sim- 
ulated spectra by iFechner et alj ([2006) for absorption lines 
in the redshift range 2.30 ^ z 2.75, who argued statistical 
comparisons must be restricted to systems with to,hi < 0.1. 

The distribution of column density ratios 77 = 
Nuaii/Nm in Fig. 1181 shows a similar decreasing trend to- 
wards increasing TVhi (top right panel). A broad distribution 
of rj values for all the systems is found (lower right panel). 
Restricting the systems to those with 0.01 < to, hi < 0.1, 
however, results in a far sharper distribution centred at 
77 ~ 1 30. The distributions are very similar to those mea- 
sured l|Fechner et al. 2006), with the distribution for sys- 
tems restricted to 0.01 < to, hi < 0.1 peaking at 77 ,$ 100 
and somewhat broader than the distribution found from the 
Enzo RT simulation. The peak in the distribution from the 
simulation is fixed by the normalisations imposed on the 
effective Hi and Hen optical depths. The agreement in the 
widths of the distribution, however, is a product of radiative 
transfer within the simulation. While the measured mean 
and breadth of the distribution may be reproduced by the 
QSO luminosity function, taking into account the Poisson 
fluctuations in QSO numbers within a Hen-i onizing photo n 
attenuation volume in a homogeneous 1GM (Meiksin 2009), 
the results found here demonstrate that a large contribution 
to the breadth in the distribution is expected to arise from 
variations in the Hen ionizing background resulting from 
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Figure 18. Comparison of Enzo RT simulation results for Hi 
and Hell line centre optical depths and column densities for fit 
absorption featres at z = 2.5. (Top left panel) Comparison of line 
centre optical depths for Hi and Hen absorption features. The 
solid line shows the relation tq jjell = To HI- (Bottom left panel) 
Hen to Hi line centre optical depth ratio tisHl line centre optical 
depth. (Top right panel) Hen to Hi column density ratio r\ vs 
Hi column density. (Bottom right panel) Distribution functions 
of r) for all lines (solid black histogram), and lines restricted to 
0.01 < tc^hi < 0.1 (dashed black histogr am). Also shown are the 
distribution functions for HS 1700+6416 llFechner et alj|2006t) for 
absorption features in the redshift range 2.30 ^ z ^ 2.75 for all 
the lines (solid red histogram) and for lines restricted to 0.01 < 
T 0,HI < 0.1 (blue dashed histogram). All distribution functions 
are normalised to unit area. 



the effect of inhomogeneities in the density field of the IGM 
on the radiation transfer of the impinging radiation field. 

The peak and FWHM in the distribution of r\ at 
z = 2.5 from the simulation correspond to an effective 
^ = Thi/ThcH — 3001}qo- For a metagalactic spectral shape 
/„ ~ u~ aMG , the ratio would be r H i/r H cii ~ 4 1+QMG . The 
match to the effective corresponds to effective spectral 
indices of 2.9 < ckmg < 3.5, much larger than the as- 
sumed source spectral index 0.5. The ratio corresponds to 
the expected softening of the inci dent spectrum due to ra- 
diative transfer through the IGM (|Madau fc Meiksin|[l994l ; 
lHaardt fc Madaul 1 1996h . The required rescaling of ^ im- 
posed on the simulations to match the measured effective Hi 
and Hen optical depths corresponds to a pre-rescaled am- 
bient radiation field with an average omg — —0.25 between 
the Hi and Hen Lyman edges arising exclusively from the 
radiative transfer of the input source spectrum through the 
simulation volume. This corresponds to the hardening effect 
of the IGM on the input radiation field. The modification 
of $ within the simulation produced by radiative transfer is 
found to be highly inhomogeneous, with typically 40 — 60 per 
cent, variations produced in the Hen ionization rate during 
the passage of the He in ionization front. The larger val- 



Figure 19. Comparison of Enzo RT simulation results for Hi 
and Hen Dopplcr parameters for fit absorption featres at Z = 2.5 
with 0.01 < tq hi < 0.1. (Top left panel) Comparison of feui an d 
^Hell absorption features. The solid line shows the relation for 
pure velocity broadening; the dashed line shows the relation for 
pure thermal broadening. (Bottom left panel) The distribution 
function of ^Hcll/frm- (Top right panel) Doppler parameter ratio 
to Hi column density. (Bottom right panel) Doppler parameter 
ratio vs Hen to Hi column density ratio r\. 

ues of j/ are comparable to those inferred from observations 
JZheng et al.ll2004 iReimers et "al]|2005l ; iFechner et al1l2006l : 
IShull et al.ll2010l ). 

The relation between the Hi and Hen Doppler param- 
eters for the absorption systems at z = 2.5 is shown in 
Fig. 1191 The Hi absorption lines are restricted to those with 
0.01 < ro, hi < 0.1. In the limit of pure thermal broad- 
ening, the Hen to Hi ion mass ratio would result in Hen 
systems with £ = &hcIi/&hi — 1/2- In the limit of velocity 
broadened lines, £ = 1. The upper left panels show that the 
Doppler parameters lie predominantly between these two 
limiting cases. Many of the points lie outside the limiting 
cases, a consequence of fitting lines to blended features as 
well as to saturated Hen features. The wide scatter suggests 
the Doppler parameter ratio £ is not a precise indicator of 
the relative contributions to line broadening due to ther- 
mal motions and internal flows. Similar results were found 
for selected unblended absorption systems in the spectrum 
of HE 2347-4342 in the redshift range 2.7 < z < 2.8, al- 
though the rat ios tended more to wards velocity broadening 
than thermal (|Zheng et al.ll2004h . No strong correlation is 
found between £ and JVhi (top right panel), however the ra- 
tio tends to increase with increasing r\ (bottom right panel), 
although this may be an effect of line deblending. 

The ratios r iiLy/ 3/r;, Lyci and r;, Ly7 /r iiLya for Hen from 
the Enzo simulation after reionization completes at z < 3 are 
shown in the upper panel of Fig. 1201 renormalizing the spec- 
tra by n^ya to the values indicated. In the region of over- 
lap, these values agree with the simulation results quoted 
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Figure 20. (Top panel) The ratios of the effective optical depths 
T l,Lyp/ T l,Lya (solid line) and Ti^ yl /ri^ ya (dashed line) for Hen 
obtained from the Enzo simulation with radiative transfer, as a 
function of Ti ^ ya . (The result is insensitive to redshift after reion- 
ization completes at z < 3.) Also shown is the predicted ratio 
T l,Lyy/ T l,hya (dotted line) based on a line-blanketing model (see 
text). (Bottom panel) The characteristic Hen t£ optical depth 
inferred from the values of Tj i^g/Tj Lya from the Enzo simulation 
with radiative transfer. 

in ISvphers et al.l ()201ll ). The inferred effective line centre 
optical depths T^ Lya using Eq. lpTjl are shown in the lower 
panel. The resulting predicted values for n^ yl /Ti^ ya from 
Eq. (|ll[) are shown by the dotted line in the upper panel. 
The values agree well with those measured directly from the 
simulation, indicating the absorption is well described by 
the line model, although the small deviation for large op- 
tical depths suggests the simple one-parameter model only 
approximately describes the underlying distribution at these 
large optical depths. 

For comparison, the value measured in the QSO 
HE 2347-4342 (z cm ~ 2.89), with r^ a ~ 2.393 ± 0.015 
at (z ) ~ 2.8, is T^Jrfft ~ 0.31 ± 0.04, using the data 
from lSyphers et al.l 1)20111 ). This agrees well with the value 
of T i.Lyp/ r i,Lya — 0.36 predicted from the Enzo simulation 
with radiative transfer, as shown in Fig. 1201 



6 COMPARISON WITH APPROXIMATE 
METHODS 

6.1 Approximate methods 

Two approximate simulation methods have been used in 
the past to characterise the effects of Hen reionization 
on the IGM. One method uses iV-body code results to 
represent the fluid by r eseating the dark m atter density 
<|Tittlev fc Meiksinll2007l ; iMcQuinn et~ai1l2009l ) ; it thus does 
not account for the hydrodynamical response of the gas 
to heating due to radiative transfer. A similar method 



post-processes hydrodynamical simulations using a radia- 
tive transfer module, but also does not include the hy- 
dro dynamical response of th e gas to the radiative trans- 
fer (|Bolton fc Haehneltil2007l ) . An alternative method artifi- 
cially boosts the Hen heating rates in optically thin gravity- 
hydrodynamics reionization simulations to mimic the boost 
in temperatures resulting from radiative transfer effects 
l|Brvan fc Machacekll2000l ). 

To estimate the importance of hydrodynamical feed- 
back effects on the observed properties of the Lya forest, 
we compare the results of the hydrodynamics simulations 
including radiative transfer with the results of the A ^-body 
radiative transfer code PMRT l|Tittlev fc Meiksin|[2007l ) . PMRT 
assumes the baryons trace the dark matter. The code is thus 
not able to cope with sudden increases in the pressure forces 
as the gas is heated and driven out of shallow potential wells 
or heated by shocks in collapsing structures. 

Past comparisons with full hydrodynamics computa- 
tions without radiative transfer show that the method re- 
covers the distribution of Doppler parameters to high ac- 
curacy, although the median Dop pler parameter may b e 
slightly broadened by 1-2 km s" 1 l|Meiksin fc WhitefeOOll ). 
A pseudo-hydrodynamics scheme with radiative transfer has 
the advantages over a fully hydrodynamical one of ease of 
implementation and a much reduced memory requirement, 
permitting simulations in larger boxes at the high spatial 
resolution required to resolve the structures in the IGM that 
give rise to the full range of measurable absorption systems. 
Here we examine how successful PMRT simulations are at cap- 
turing the essential features of the Lya forest predicted by 
full hydrodynamics computations including radiative trans- 
fer. 

The PMRT runs used 512 3 gravitating particles, with the 
identical initial particle distributions in position and velocity 
as for the Enzo runs, and the identical cosmological param- 
eters. It used the identical radiative transfer module as for 
the Enzo runs, using 256 3 cell s to represent the fluid. Follow- 
ing |Tittle^^Miiksi3 (|2007l ). all cells with fewer than 100 
particles were convolved over a radius of 1.5 cells to match 
the gas distribution from Enzo. 

We also implement the alternative method of artificially 
boosting the Hen heating rate in a hydrodynamical simula- 
tion in the optically thin limit to make up for the shortfall 
in the additional heating due to radiative transfer effects. A 
boost by a factor of 2 was found to provide sufficient heat to 
broade n the absorption lines clos e to the measured values at 
z ~ 3 (|Brvan fc Mac hacek 200(3). We recompute the reion- 
ization simulation using Enzo with the Hen heating rate per 
Hen ion boosted to G Ho ii = 2Gh c ii, and the radiative trans- 
fer switched off. 

6.2 Physical properties of IGM 

Maps of the IGM temperature are shown in Fig. [2T] at 
z = 3.3, shortly before Hen reionization completes, and in 
Fig. [22] at z = 2, well after completion. At z = 3.3, the 
PMRT run matches most closely to the Enzo with RT sim- 
ulation. There are, however, some conspicuous differences. 
The temperature structures tend to be sharper in the PMRT 
simulation, as it lacks both shock heating and the pressure 
forces that smooth the gas. The effect of shadowing by dense 
clumps is also more prevalent in the PMRT simulation. This is 
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Figure 21. Temperature map at z = 3.3, shortly before Hen 
reionization completes. The ionization front sweeps across the 
boxes from top to bottom. The maps are shown for Enzo with 
radiative transfer (upper left), in the optically thin limit (upper 
right), in the optically thin limit with the Hell heating rate dou- 
bled (lower right), and for the pseudo- hydrodynamics code PMRT 
(lower left). The box side is 25h~ 1 Mpc (comoving). 



Figure 22. Temperature map at z = 2.0, well after Hell reion- 
ization completes. The ionization front swept across the boxes 
from top to bottom. The maps are shown for Enzo with radiative 
transfer (upper left), in the optically thin limit (upper right), in 
the optically thin limit with the Hen heating rate doubled (lower 
right), and for the pseudo-hydrodynamics code PMRT (lower left). 
The box side is 25/i _1 Mpc (comoving). 



a consequence of the inability of the dense regions to expand 
due to pressure forces as the ionization front sweeps across 
and boosts the gas temperature and pressure, as discussed 
in Sec. 14.21 The Enzo simulation with the boosted Hen heat- 
ing rate recovers the high temperatures in the filaments, but 
the temperature in the underdense regions, although warmer 
than the optically thin simulation, still lies well below the 
Enzo simulation with radiative transfer. 

By z — 2.0, the PMRT simulation shows significant dis- 
crepancies in the gas temperature within the filaments com- 
pared with the Enzo RT simulation. In contrast to the hydro- 
dynamics simulation which accounts for the hydrodynami- 
cal response of the gas to the extra heating following Hen 
reionization, allowing it to expand and cool, the PMRT simu- 
lation keeps the gas too hot. In underdense regions, the PMRT 
simulation over-estimates the amount of adiabatic cooling, 
resulting in the gas becoming too cold. The Enzo simula- 
tion with the boosted Hen heating rate is better able to 
reproduce the temperatures in the underdense regions, but 
maintains the gas in the filaments at too high temperatures 
compared with the Enzo RT simulation. 

These effects may be quantified using the temperature- 
density distributions extracted from the simulations. The 
temperature-density distribution obtained with PMRT closely 
resembles that found using Enzo with radiative transfer at 
z > 3.0, as shown in Fig. 1231 with reasonably good agree- 
ment down to z — 2. The distributions overlap very closely, 
although the hydrodynamics computation shows a some- 
what tighter and flatter distribution in temperatures, espe- 
cially towards z = 2. While the temperatures in the optically 
thin simulation with boosted Hen heating rate agree well 
with the RT simulations for p/{p) > 3 over 3.0 < z < 3.5, 
the temperatures are appreciably lower in structures with 
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Figure 23. Temperature-density distributions at z = 4.0, 3.5, 3.0 
and 2.0 following Hell reionization, for computation using Enzo 
with radiative transfer, Enzo in the optically thin limit, Enzo in 
the optically thin limit but with the Hell heating rate doubled, 
and the pseudo-hydrodynamics RT code PMRT. Note the smaller 
temperature scale at z = 4.0, prior to Hell reionization. The 
distributions are normalised to unity. The contour levels are at 
probability density levels (per d log 10 T-d log 10 p/{p)) 0.1, 0.32, 
1.0, 3.2, 10, 32, ... . 




Figure 24. Temperature distributions at z = 4.1, 3.5, 3.0 and 2.0 
for computation using Enzo with radiative transfer (solid lines), 
Enzo in the optically thin limit (dashed lines), Enzo in the op- 
tically thin limit but with the Hen heating rate doubled (dot- 
dashed lines), and the pseudo- hydrodynamics RT code PMRT (dot- 
ted lines). Note the smaller temperature scale at z = 4.1, prior to 
Hen reionization. 

p/{p) < 1- At = 2, good agreement is found between the 
Enzo RT and boosted Hen heating rate simulations for un- 
derdense structures, but the latter produces too high tem- 
peratures for overdense structures. The distribution from 
the optically thin Enzo simulation agrees well with the oth- 
ers at z = 4, but once Hen reionization commences, the 
temperatures are systematically too low. 

The temperature distributions of the Enzo with radia- 
tive transfer and the PMRT runs are very similar, as shown in 
Fig. 1241 which may be expected from the overall agreement 
in the density-temperature distrubtions. By z — 2.0, how- 
ever, the PMRT simulation systematically underestimates the 
temperatures by 3 — 5 x 10 3 K. Most of the volume is too cold 
in the Enzo simulation with the boosted Hen heating rate 
for 3.0 < z < 3.5, although a tail of warmer temperatures 
results from overdense material, as shown in Fig. [23] By 
z — 2, these temperatures extend the distribution towards 
too high values compared with the RT simulations. 

6.3 Spectral signatures 

To assess the observational impact of the physical differ- 
ences arising from the hydrodynamical response of the gas 
compared with the pseudo-hydrodynamics simulation on the 
Lycv forest, sample spectra are drawn from the Enzo and 
PMRT simulations, and the absorption features fit. As is clear 
from Figs. [8] and [9] reionization induces substantial peculiar 
velocities to the gas that a pseudo-hydrodynamics scheme 
like PMRT is unable to reproduce. The differences are small 
compared with the peculiar velocities present due to large 
scale flows, as is evident by a comparison with the full pe- 
culiar velocity distribution shown by the rightmost curve in 
Fig. [9] The comparison between the Enzo run with radiative 
transfer and the run without reionization, however, suggests 
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Figure 25. Evolution in the median Doppler parameter for the 
computation using Enzo with radiative transfer (solid lines), Enzo 
in the optically thin limit (dashed lines), Enzo in the optically thin 
limit with the Hell heating rate doubled (dot-dashed lines), and 
the pseudo- hydrodynamics code PMRT (dotted lines). (Top panel) 
Results for lines optically thin at line centre. (Bottom panel) Re- 
sults for lines optically thick at line centre. 



differences as great as 10 km s _1 are produced, with differ- 
ences of 1 — 2 km s _1 typical. These should be reflected in 
the line widths, although the differences in gas temperature 
will affect the line widths as well. For a pre-Hen ionization 
temperature of 5000 K, a temperature boost to 20000 K will 
double the thermal line width, increasing it to 18 km s _1 . 
The additional broadening due to the increase in peculiar 
velocity will augment the increase due to thermal broaden- 
ing by at most only a few kilometres per second. In com- 
parison with the optically thin reionization simulation, the 
additional peculiar velocities induced when radiative trans- 
fer is included add negligibly to the broadening. 

The evolution of the resulting median Doppler parame- 
ters is shown in Fig. [25] The values generally agree well for 
the absorption lines optically thick at the line centre (bot- 
tom panel), although with deviations of ~ ±4 km s _1 from 
the Enzo RT simulation. The optically thin Enzo simulations 
agree best with the RT run at z > 3, with the simulations 
with the boosted and non-boosted Hen heating rates brack- 
eting the RT simulation, but producing too broad features 
by z < 3, where the PMRT simulation gives near perfect agree- 
ment. 

By contrast, the PMRT simulation recovers the median 
Doppler parameters for the absorption lines optically thin at 
line centre for z > 3. These weaker lines arise from gas too 
rarefied to maintain atomic thermal balance with the pho- 
toionization heating rate, so that including the temperature 
heating due to radiative transfer is required. The good agree- 
ment with the Enzo RT simulation results shows that hydro- 
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dynamical feedback constributes little additional broadening 
to the lines. The optically thin Enzo simulation produces 
lines that are ~ 5 km s" 1 too narrow at 3 < z < 3.5. Boost- 
ing the Hen heating rate recovers the median Doppler pa- 
rameter well for 3 < z < 3.5, but overheats the gas at lower 
redshifts. At z < 3, adiabatic cooling in the PMRT simulation 
produces features that are too narrow compared with the 
Enzo RT simulation. None of the approximate methods pro- 
vides an adequate description of the absorption line widths 
over the full redshift range of interest. 



7 CONCLUSIONS 

A radiative transfer m odule based on the method of 
iTittlev fc Meiksinl l|2007l) was incorporated into the gravity- 
hydrodynamics code Enzo to investigate the impact of Hen 
reionization on the structure of the IGM. The RT module 
is specifically designed to account for the heating by the 
hardening of the radiation field on passing through inho- 
mogeneities in the density field, crucial for capturing the 
heating due to hard sources like QSOs. The combined code 
is able to capture the hydrodynamical feedback following 
the reionization of both hydrogen and helium. 

A simulation of a plane-wave ionization front sweeping 
across a volume 25ft _1 Mpc (comoving) on a side, repre- 
senting a characteristic /-front produced by photoionizing 
sources, was performed. The source spectrum was initially 
a starburst, transforming gradually into a QSO spectrum 
f v ~ u~ ' 5 over the redshift interval 3 < z < 4 to repre- 
sent the onset of Hen reionization by QSO sources. Sim- 
ulations in the optically thin limit, one without and one 
with the Hen heating rate doubled, and an JV-body based 
pseudo-hydrodynamics simulation with radiative transfer 
using PMRT, are performed with identical initial conditions to 
the Enzo radiative transfer simulation in order to examine 
how well results from these simpler approximate methods 
compare. 

Substantial heating is introduced by Hen reionization, 
with temperatures exceeding those following Hen reioniza- 
tion in the optically thin limit by typically 10 — 20 x 10 s K. 
High temperatures, in excess of 35 x 10 3 K, are pro- 
duced in the "shadows" of overdense structures as the Hen 
is photoionized by a radiation field hardened on passing 
through the overdense structures. While peculiar velocities 
are boosted by 0.1 — 10 km s _1 by z = 3 over a simulation 
with no Hen reionization, the differences compared with a 
simulation with Hen reionization in the optically thin limit 
are typically smaller than 1km s _1 . The expansion of the 
gas following Hen reionization, however, plays an impor- 
tant role in curtailing the hardening of the radiation field 
by reducing the gas density in small haloes. As a result, the 
temperature-density relation is somewhat tighter and flat- 
ter compared with the corresponding PMRT simulation, for 
which no hydrodynamical response is accounted for. In rare 
regions, the temperature in dense regions when radiative 
transfer is included is lower than the corresponding Enzo 
simulation in the optically thin limit, a consequence of com- 
plete shadowing of Hen- ionizing photons by dense clumps 
in the radiative transfer computation. 

The gas flows induced by Hen reionization alter the gas 
density field on large scales (several comoving megaparsecs 



in extent), with density reductions by as much as 10-20 per 
cent, compared with a simulation without Hen reionization. 
In much smaller regions, the gas is pushed into structures 
with density increases of up to 50 per cent. Compared with 
a simulation with Hen reionized in the optically thin limit, 
the density changes are much reduced, with large scale re- 
ductions of ~ 5 per cent., and density enhancements of a 
few to several per cent, in complex density regions. 

Because the dark matter is coupled gravitationally to 
the gas, displacements in the gas will produce displacements 
in the dark matter as well. Large scale coherent reductions 
in the dark matter density of up to 0.5 per cent, are found at 
z = 3, increasing to nearly 1 per cent, by z = 2. Small dense 
regions show enhancements in density by 0.2 per cent, at z = 
3 to 1 per cent, at z = 2. The dark matter and gas density 
displacements may produce distortions of a few per cent, in 
the matter power spectrum on comoving wavenumber scales 
of k > 0.5/iMpc -1 . 

The reionization of Hen has a substantial impact on the 
Hi Lya forest, as would be measured in the spectra of back- 
ground QSOs. Synthetic spectra are drawn from the simu- 
lations and normalised according to measured values of the 
mean Hi and Hen transmission values. The most conspic- 
uous effect is on the broadening of the absorption features. 
As the Doppler parameter distribution may be measured 
to high accuracy in a given QSO spectrum, with a median 
Doppler parameter determined to a precision of ~ 1 km s _1 
or better, the line widths are a potentially powerful means 
of constraining the Hen reionization process and the nature 
of the ionizing sources. (Alternative measures of the line 
shapes, such as discrete wavelet coefficients or the flux cur- 
vature, may similarly be sensitive probes of the reionization 
history.) The heating induced by radiative transfer, however, 
will affect the other flux-related statistics as well. 

Compared with Hen reionization in the optically thin 
limit, the boost in temperatures when radiative transfer is 
included produces a larger fraction of weak Hi absorption 
features with optical depths < 0.1, and a smaller fraction of 
moderate optical depth systems with optical depths 0.1 — 0.6. 
This may be accounted for by the increase in the line widths 
and reduction in the neutral fraction at a given gas den- 
sity due to the higher temperatures, as both tend to reduce 
the line centre optical depth of the absorption features. The 
modifications to the optical depths result in a shift in the 
pixel flux cumulative distribution, with a maximum differ- 
ence between the cumulative distributions from the Enzo 
runs with and without radiative transfer of ~ 0.1 at z = 3, 
readily detectable in measured spectra. 

The Hi Lya forest spectra are found to be fully resolv- 
able into absorption features. The resulting Hi column den- 
sity distribution of the fit lines is somewhat broader when 
radiative transfer is taken into account compared with the 
optically thin simulation. Fewer narrow lines are found in 
the simulation with radiative transfer than without. The me- 
dian Doppler parameter evolves little between 2.5 < z < 3.5, 
with femed — 35 km s _1 , but decreases to 6 me d 25 km s _1 
by 2 = 2. Dividing the absorption features into those opti- 
cally thin and optically thick at line centre reveals a sharp 
difference in the evolution of the features. Those optically 
thick evolve slowly, increasing from fo mc d — 27 km s _1 when 
Hen reionization begins at z = 4.0 to 6 mo d — 33 km s _1 at 
z = 3.5, then gradually diminishing to & mc d — 19 km s" 1 at 
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z = 2.0. By contrast, the median Doppler parameter of the 
optically thin lines increases rapidly from 6 me d — 21 km s _1 
at z = 4.0 to 6 mcd ~ 44 - 45 km s" 1 at 3.0 < z < 3.5, 
then diminishes slightly to 6 me d — 39 km s -1 at z — 2.5 
and rapidly to 6 mc d — 27 km s _1 by z = 2.0. The me- 
dian Doppler widths of the optically thick systems agree 
well with those measured. By contrast, the median values 
of the optically thin systems at 3 < z < 3.5 substantially 
exceed those measured, suggesting either Hen reionization 
was imposed too late in the simulation, having just com- 
pleted at 2 ~ 3, or that the source spectrum is too hard. A 
previous PHRT simulation with He n reionization initiated at 
z = 5.0 l|Tittlev fc Meiksmll2007l ) produced a substantially 
lower peak temperature of 15 x 10 3 K by z = 3.0, cooler 
than found here by 10 — 20 x 10 3 K. Cooler temperatures 
may thus be achieved if Hen reionization were substantially 
completed prior to z = 3, although this may be at tension 
with recent claims that Hen reionization extended to z < 3. 
The PMRT simulations also found cooler temperatures arose 
when a single power-law spectrum reionized both hydrogen 
and helium, so that the temperatures are sensitive to the 
amount of Hen reionization during the hydrogen reioniza- 
tion epoch. Clearly simulations allowing for a more complex 
Hen reionization history involving a range of source histories 
and spectra are required to fully match the data. 

The PMRT simulation matches the median Doppler width 
of the optically thick absorbers in the Enzo radiative trans- 
fer simulation to an accuracy of 2 km s _1 at z > 4 and 
better than 0.5 km s _1 at 2 < z < 2.5, but underpredicts 
the widths by as much as 3.5 km s~ at intermediate red- 
shifts. For the absorbers optically thin at line centre, the 
agreement is to better than 2 km s _1 for 3.0 < z < 5, but 
underpredicts the widths by ~ 5 km s -1 by z = 2.0. 

Doubling the Hen heating rate in an optically thin Enzo 
simulation produces a close match to the median Doppler 
parameter of the optically thin absorbers from the Enzo sim- 
ulation with radiative transfer over 3.0 < z < 3.5, agreeing 
to within 0.6 km s _1 . But by z = 2.0, the median Doppler 
parameter is overpredicted by ~ 15 km s _1 . The Doppler 
parameters of the absorption systems optically thick at line 
centre are overpredicted by 2-6 km s _1 over 2.0 < z < 3.5. 
The approximate methods are thus unable to reproduce the 
correct amount of broadening to the accuracy with which it 
is measured over the full range of redshifts for which the for- 
est may be observed in high resolution, high signal-to-noise 
ratio spectra. 

Just as absorption features may fully account for the 
pixel flux distribution in Lya, a line blanketing model ac- 
counts well for the evolution in the mean transparency of 
the IGM in higher order Lyman series lines. Effective op- 
tical depth ratios for the Hi of TLys/^Lya ~ 0.36 and 
r Ly7/ r Lya — 0-21 are predicted at z = 3.0 from the Enzo 
radiative transfer simulation. 

Radiative transfer through density inhomogeneities in 
the IGM results in a wide spread of Hen Lya optical depths 
compared with Hi even after Hen reionization has com- 
pleted in the simulation at z = 3. A broad peak centred 
at thcIi/thi — 105 is found at z = 3.0, diminishing to 80 at 
z = 2.5 and 75 at z = 2.0. 

As for the Hi, the pixel flux distribution for Hen Lya 
may be fully accounted for by absorption lines for 2.0 < 
z < 2.5. At z = 3.0, however, the mean transmission is too 



low for meaningful absorption line fitting to be performed. 
Matching to the corresponding Hi absorption features re- 
sults in a broad distribution of the column density ratio 
r\ — -/VhcIi/Ahi, with a peak of r\ ~ 130 at z = 2.5, well 
matching, though somewhat narrower than, the distribution 
measured in the ultraviolet spectra of high redshifts QSOs. 
While the value for 77 at the distribution peak results from 
the rescaling of the relative Hi and Hen ionization rates re- 
quired to match the measured effective Hi and Hen Lya 
optical depths, the distribution width is produced by the 
simulation and suggests that much of the observed spread 
may arise naturally from the radiative transfer through an 
inhomogeneous IGM. 

Comparison between the column density and optical 
depth ratios indicates a typical Doppler parameter ratio 
of £ = &Hcii/&Hi — 0.4, close to the thermally broad- 
ened limit of f = 0.5. A comparison based on well de- 
blended, unsaturated systems with line centre Hi optical 
depths 0.01 < ro < 0.1 shows a range in Doppler parameter 
ratios peaking between 0.5 < £ < 1.0, the latter correspond- 
ing to the limit of velocity-broadened features. 

Line blanketing accounts well for the higher order Hen 
Lyman transmission through the IGM as well. The simula- 
tions predict effective optical depth ratios of i"Ly(3/ T Lyn — 
0.36 at z = 2.5, in agreement with the measured value, and 
T Ly7/ r Lya — 0.20 (currently unmeasured). 

A range of models, allowing for multiple QSO sources 
turning on over a range of redshifts with a variety of lumi- 
nosities, spectral shapes and lifetimes, would be required to 
accurately match the Hi and Hen absorption data as mea- 
sured in the spectra of high redshift background QSOs. The 
simulations presented here demonstrate that Hen reioniza- 
tion simulations including radiative transfer are required to 
fully exploit the growing body of high quality IGM data at 
both optical and ultraviolet wavelengths. 
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